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Abstract 

All regular and singular cosmological perturbations in a radiation dominated Einstein- 
de Sitter Universe with collisionless particles can be found by a generalized power 
■^^ ■ series ansatz. 

Although the study of cosmological perturbations has a long history iQ , exact solutions 
in the linear regime are known for a few perfect fluid models only (see, e.g., pi). Within a 
more fundamental description of matter no such solutions have been found until recently 
[pi H, p|. In its early epochs the Universe contains various collisionless, massless particles. 
Examples are background gravitons generated during the Planck epoch and neutrinos 
after their decoupling from the electrons. 

Cosmological perturbations for collisionless matter have been analyzed numerically 

(e.g., m) and analytically in the short- and long- wavelength approximations M. Two 

^ ■ kinds of solutions show up: Regular solutions are specified by bounded metric pertur- 

t^^ I bations at the Big Bang singularity, whereas for singular solutions the metric blows up, 

w^ ■ changing the type of the initial singularity. The latter show superhorizon oscillations 

as discussed by Zakharov and Vishniac [g|. The general solution is a combination 

V-v ■ of both types. Only the regular part of the solution depends on the distribution of the 

T^lj- I collisionless particles at the singularity. 

0^ ■ We obtain exact solutions for scalar perturbations by a generalized power series 

fj I ansatz. The matter under investigation is a mixture of collisionless particles and a per- 

^^i feet fluid in a radiation dominated Einstein-de Sitter background. The corresponding 

discussion for the vector and tensor perturbations can be found in [§[ g[. The cqua- 
bX)i tions of motion are formulated gauge invariantly. The matter perturbations can either 

be obtained from the Boltzmann equation ||7|, |lO| or from the flnite-temperature graviton 
self-energy Q. Reference Q compares the two approaches. 

We use Bardeen's gauge invariant metric potentials |1^. This is useful since the 
S i perturbation of the energy-momentum tensor 



X 



is gauge invariant from the deflnition of the effective action F*^ (for background-covariant 
gauges). The scalar metric potentials are denoted by $ and H. They are related to the 
gauge invariant density contrast e™ (which is the density contrast 6 in the local rest-frame 
of matter) and the anisotropic pressure ttt by the Einstein equations: 

Y'^ = e™, (2) 

X^n = TTT . (3) 

The variable x is the product of the conformal time and the (fixed) wavenumber of the 
perturbation. 
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The high-temperature hmit of the graviton self-energy (the second varation of T^ in 
Eq. (n])) for temperatures below the Planck scale and the kinetic theory yield the same 
equations of motion [|| for the metric potentials. They form a closed set of equations 
allowing isentropic perturbations in the perfect fluid component only. With a being 
the ratio of the energy density in coUisionless matter and the total energy density, and 
^ = ^cii + ^PF, the trace of the Einstein equations reads: 
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The 7„'s are related to the n-th momenta of the coUisionless particles' initial distribution 
function [g. The evolution starts at the Big Bang (x — 0). 
To solve the equations (|4|) - (g) we make the ansatz 



F{x) = CiF„g(x) + C2x"F,i„g(a;) 
where Ci and C2 are arbitrary constants and 



(7) 



rcg,sing:^ 
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The singular solutions do not depend on the initial conditions 7„ . They can be added to 
any regular solution specified by the initial conditions. The exponent 



cr± 



± 



5 -32a 



20 



(9) 



takes complex values for a > acrit — 5/32. Thus superhorizon size oscillations occur. 
They can not be taken seriously up to the initial singularity, because nonlinear effects of 
the Einstein equation will be important. But they are necessary to match a mixmaster 
scenario [|l2| or quantum generated perturbations from a preceding inflationary epoch 
[ [13| . A strictly Friedmannian singularity rules out the singular solutions. 

In the Figs. 1 and 2 we consider a model with three massless neutrino species and a 
perfect fluid consisting of relativistic electrons and photons in equilibrium. For such a 
model a- 1/2 |l|. 

Fig. 1 shows the regular solutions for 11(0) = and $(0) — 1, $'(0) = 0, i.e., a perfect 
fluid at the origin. The density contrast grows on superhorizon scales and oscillates after 
horizon crossing [x/tt = 1). On subhorizont scales the perfect fluid component has con- 
stant amplitude, whereas the coUisionless component decays (^ I/2;) due to directional 
dispersion. 
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Figure 1: Regular scalar perturbations in a two component universe (a — 1/2). The 
various lines show regular solutions for \5cu\ (full line), \bpF\ (dotted line), and \'kt\ 
(dashed line). The initial conditions are specified by 70 == 5/64, 72 — 427o, and 74 — 497o. 
All other 7„'s vanish. 




Figure 2: Singular solutions corresponding to oj^. Same variables as in Fig. 1. Notice 
the logarithmic scale in x — the solutions are essentially singular as a; ^ 0. The second 
solution o- differs by a phase. 



In Fig. 2 the singular solutions are plotted. The perfect fluid has the same behaviour 
as in Fig. 1, but the coUisionless component shows superhorizon oscillations as mentioned 
above. These oscillations are of geometrical, not of causal, origin. Lukash et al. showed 
that coUisionless matter can isotropize an initial anisotropy p5[ . Reversing time these 
solutions reflect the instability of an isotropic collapse. The sound velocity inside the 
horizon equals to the speed of light for all solutions. 
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